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Abstract. We propose a time-space discretization of a general notion of quasistatic 
growth of brittle fractures in elastic bodies proposed in by G. Dal Maso, G.A. 
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1. Introduction 

The aim of this paper is to provide a discontinuous finite element approximation of a model of 
quasistatic growth of brittle fractures in finite elasticity recently proposed in ^21 by Dal Maso, 
Francfort, and Toader in the framework of the variational theory of crack propagation proposed 
by Francfort and Marigo in 15 . This theory is inspired to Griffith's criterion and determines the 
crack path through a competition between bulk and surface energies. 

In the case of linearized elasticity, a first precise mathematical formulation of the model jl5| 
has been given by Dal Maso and Toader jl2) : they treat the case of anti-planar shear in dimension 
two assuming that the fractures are compact sets with a finite number of connected components. 
This analysis has been extended to the case of plane elasticity by ChamboUe in P. Francfort 
and Larsen ^1], using the framework of SBV functions (see Section [21), proved the existence of 
a quasistatic growth of brittle fractures in the case of anti-planar shear in any dimension N > 2 
and without assumptions on the structure of the fractures which are dealt with the set of jumps 
of the displacements. Approximation results for the quasistatic evolution of |14| has been given in 
|16| and in |17| and provide a theoretical basis to the munerical study of the model given in [S]. 

The quasistatic crack growth proposed by Dal Maso, Francfort, and Toader in consider 
the case of finite elasticity, and takes into account possible volume and traction forces applied to 
the elastic body. In order to describe the result of 13 (a complete description is given in Section 
O, let us assume that the elastic body has a reference configuration given by C open, 
bounded and with Lipschitz boundary. Let QdJI C dfl be open in the relative topology, and let 
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dN^l :— d^\dj:)il. Let C fl, and let ds^ ^ d^fl be such that fi^ <^ds^ ~ 0- is the brittle 
part of f2, and i9sr2 is the part of the boundary where traction forces are supposed to act. A crack 
is given by any rectifiable set in with finite {N — 1) Hausdorff measure. Given a boundary 
deformation g on dn^ and a crack F, the family of all admissible deformation of $1 is given by 
the set AD{g,r) of all function u € GSBV{n;R'^) (see SectionEJ such that S{u) C T and m = 5 
on do^ \ r. Here S{u) denotes the set of jumps of u, and the equality it = 5 is intended in the 
sense of traces. Requiring u ^ g only on \ T means that the deformation is assumed not to 
be transmitted through the fracture. The bulk energy considered in '13' is of the form 

W{x,Vu{x))dx, 

where W{x,S,) is quasiconvex in and satisfies suitable regularity and growth assumptions (see 
(E3I) and lESIl)- Moreover the time dependent body and traction forces are supposed to be 
conservative with work given by 

-/ F{t,x,u{x))dx - I G{t,x,u{x))dn^~^{x), 
Jn\r J dsn 

where F and G satisfy suitable regularity and growth conditions (see Section EJ. Finally the work 
made to produce the crack F is given by 

S'iT) := J kix,iy{x))dn^-\x), 

where I'ix) is the normal to F at x, and k{x, v) satisfies standard hypotheses which guarantee 
lower semicontinuity (see Section OJ. Clearly, W,F,G and k depend on the material. Let us set 

£''\t){u):= I W{x,Vu{x))dx- [ F{t,x,u{x)) dx - [ G{t,x,u{x))dH^-^{x), 

and 

(1.1) £{t){u,T) ■.= £'''{t){u)+£'{T). 

Given a boundary deformation g{t) with t £ [0,T] and a preexisting crack Fq, a quasistatic crack 
growth relative to g and Fq is a map {t (ii(i),F(t)) : i e [Oi^^]} such that the following 
conditions hold: 

(1) for all t e [0,r]: u(t) e AD{g(t),r{t)y, 

(2) irreversibility: Fq C F(s) C F(i) for all < s < i < T; 

(3) static equilibrium: for all t G [0, T] and for all admissible configurations (m, F) with T(t) C F 

£{t)iu{t),T{t))<£it){u,T); 

(4) nondissipativity: the time derivative of the total energy £{t){u{t),T(t)) is equal to the 
power of external forces (see H3.21|l '). 

In this paper we discretize the model using a suitable finite element method and prove its con- 
vergence to this notion of quasistatic crack growth. We restrict our analysis to a two dimensional 
setting considering only a polygonal reference configuration f2 C R^. 

The discretization of the domain fl is carried out as in jl7| employing adaptive triangulations 
introduced by M. Negri in 19 (see also 20J). Let us fix two parameters £ > and a g]0, We 
consider a regular triangulation of size e of Q, i.e. we assume that there exist two constants Ci 
and C2 so that every triangle T G contains a ball of diameter ciE and is contained in a ball of 
diameter C2e. In order to treat the boundary data, we assume also that dn^ is composed of edges 
of Rj. On each edge [x,y] of Rg we consider a point z such that z — tx+{l — t)y with t S [a, 1 — a]. 
These points are called adaptive vertices. Connecting together the adaptive vertices, we divide 
every T S Re into four triangles. We take the new triangulation T obtained after this division as 
the discretization of The family of all such triangulations will be denoted by T^.ai^)- 
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The discretization of the energy functional is obtained restricting the total energy Hl.l|) to the 
family of functions u which are afhne on the triangles of some triangulation T(m) e 7^_a(r2) and are 
allowed to jump across the edges of T(u) contained in fi^. We indicate this space by Af^a{Vl] M?). 
The boundary data is assumed to belong to the space ATe{^]^^) of continuous functions which 
are afhne on every triangle T G R^. 

Given the boundary data S W^'^{[0,T],W^^p{Vl-S?) CMtiVL-M?)) with g^it) £ AFe{^]R^) 
for all t e [0, T] {p, q are related to the growth assumptions on W, F, G) and an initial crack r° ^ 
(see SectionEl, we divide [0, 1] into subintervals [if , tf+i] of size (5 > for i = 0, . . . , Ns, and for all 
u £ ^^^(fi;R^) we indicate by S^^^\u) the edges of the triangulation T{u) contained in dn^ on 
which u ^ (?e(i)- Using a variational argument (Proposition!^^, we construct a discrete evolution 
{(wf'j,, T^'^J : i = 0,...,Ns} such that for alH = 0, . . . , TV^ we have ttf;j, € A^^^i^; R"^), 

i 

and the following unilateral minimality property holds: for all v E ■4^^(f2;R^) 
(1.2) S^'(tt){4-^ < £^'{tf)iv)+£^ ((^(z;) U^^^(*-^(«)) \rf;r^) ■ 

Notice that by construction uf^^ G AD{g^{tf), Pf'^). Moreover the definition of the discrete fracture 
ensures that Pf C Pf.J for all i < j, recovering in this discrete setting the irreversibilty of the crack 
growth given in (2) . The minimality property p.2|l is the reformulation in the finite element space 
of the equilibrium condition (3). Finally we obtain an estimate from above for £(<*)(u^'^, P*'^) 
(see Proposition l6.2(l which is a discrete version of (4). 

In order to perform the asymptotic analysis of the discrete evolution {(Wg'^,P^'^) : i — 
0, . . . ,Ns} we make the piecewise constant interpolation in time u^^^{t) = wf'^ and Pf jj(i) = Pf'_^ 
for all tf <t < tf_^j. Let us suppose that 

5e ^ g strongly in W^''^{[0,T],W^'P{n;R'^) n L^n;^^)) 

(where on VP^'^'(ri; R^) n L'^(il;IR^) we take the norm ||m|| := ||u||h'1'P(0;R2) + ||w||l9(0;R2)), and 
that Pg ^ approximate an initial crack P" in the sense of Proposition 15. II 

The main result of the paper (Theorem I7.1|l states that there exist a quasistatic evolution 
{t (w(i), r(t)) : t E [0,T]} in the sense of ^31 relative to the boundary deformation g and the 
preexisting crack P° and sequences (5„ — *■ 0, e„ ^ 0, a„ — + 0, such that setting 

u,,it):^uil^^Jt), P„(i):=Pf;;,„Jt), 

for all t e [0, T] the following facts hold: 

(a) (u„(t))ngN is weakly precompact in GSBV^{il;M.'^), and every accumulation point u{t) is 
such that u{t) e AD{g{t), P(t)), and {u{t), T{t)) satisfy the static equilibrium (2); moreover 
there exists a subsequence (Sn^, Enj., anfc)fceN of (J„, £„, a„)„gN (depending on t) such that 

{t) u{t) weakly in GSBV^{n; K^) 

(see Section !2| for a precise definition of GSBV^{Vl\'R?), and of weak convergence in this 
space); 

(b) convergence of the total energy holds, and more precisely elastic and surface energies 
converge separately, that is 

E^\t){ur,{t)) S^\t){u{t)), 5^(P„(t)) -> s'{r{t)). 

By point (a), the approximation of the deformation u(t) is available only up to a subsequence 
depending on t: this is due to the possible non uniqueness of the minimum energy deformation 
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associated to T{t). In the case E'^^{t){u) is strictly convex, it turns out that the deformation u{t) 
is uniquely determined, and we prove that fTheorem l8.1|) 

Vu„(t) -> Vu{t) strongly in LP{Q.] M'^''^), 

and 

u„{t) u{t) strongly in L'^{n;R^). 

The main difficulty to prove Theorem 17. II consists in passing to the limit in the static equilibrium 
(|1.2|) . In order to find the fracture T{t) in the limit, in Lemma [7.21 and Lemma f7. 41 we adapt to 
the context of finite elements the notion of a^- convergence of sets formulated in [TH]. This is 
the key tool to obtain the convergence of elastic and surface energies at all times t e [0, T] (while 
in this was available only at the continuity points of 7i^(r{t))). In order to infer the static 
equilibrium of T(t) from that of r„(t), we employ a generalization of the piecewise affine transfer 
of jumps 23 Proposition 5.1] (see Proposition I4.2|l . 

The paper is organized as follows. In Section[21we introduce the basic notation, and some tools 
employed throughout the paper. In Section |21 we describe the quasistatic crack growth of |13| 
precising the functional setting and the hypotheses on the elastic and surface energies involved. In 
Section^lwe introduce the finite element space, and in Section[3we prove an approximation result 
for a preexisting crack configuration. In Section we prove the existence of a discrete evolution, 
and in Section [7| we prove the main approximation result (Theorem 17. In Section |S1 we treat 
the case of strictly convex total energy. 

2. Notations and Preliminaries 

In this section we introduce the main notations and the preliminary results employed in the 
rest of the paper. 

Basic notation. We will employ the following basic notation: 

- ][{nxm -j-jj^g space of n X TO matrices; 

- Ti,^ is the one-dimensional Hausdorff measure; 

- for p G [1, +oo], II • lip denotes the usual norm; 

- if /i is a measure on and A is a Borel subset of M^, fil—A denotes the restriction of jj. 
to A, i.e. (^Lv4)(B) := ^i{B n A) for all Borel sets B C M^. 

- if A, B C M.^ ^ A C B means that A C S up to a set of Ti^-measure zero. 



SBV and GSBV spaces. Let A be an open subset of M", and let u : A — > be a measurable 
function. Given a; e A, we say that u{x) is the approximate limit of m at a;, and we write 
u{x) = aplimu(j/), if for every e > 

lim r'"/:" {{y G B,.(x) : \u{y) - u{x)\ > e}) = 0. 

Here Br{x) denotes the ball of center x and radius r. We indicate by S{u) the set of points where 
the approximate limit of u does not exist. We say that the matrix mxn Vu(a;) is the approximate 
gradient of m at a; if 

ap lim <y)~<^) ~^<^)iy-^) ^ 0. 

y^x \y-x\ 

We say that u e BV{A;M."^) if u G L^{A;R"^), and its distributional derivative Du is a vector- 
valued Radon measure on A. In this case, it turns out that S{u) is rectifiable, that is there exists a 
sequence (Afi)igN of C^-manifolds such that S{u) C |J. Ali up to a set of '^"^^-measure zero; as a 
consequence S{u) admits a normal i/^ for 7i"^^-almost every x G S{u). Moreover the approximate 
gradient Vu(a;) exists for a.e. x € A, and Vu is the density of the absolutely continuous part of 
Du. 

We say that u G SBV{A; M™) if u G BV{A; M™) and the singular part D^u of its distributional 
derivative Du is concentrated on S{u). The space SBV{A;M.'^) is called the space of R™-valued 
special functions of bounded variation. For more details, the reader is referred to |4i. We indicate 
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with SBViociA^"') the space of functions which belong to SBV{A' ,W"') for every open set A' 
with compact closure in A. 

The set GSBV{A,W") is defined as the set of functions u : A ^ such that (p{u) G 
SBVioc{A) for every (f G C^(R'") such that the support of has compact closure in R™. If 
p g]1, +oo[, we set 

GSBVP{A,R"') {u G GSBViA^M."") : Vu G Lf(A,M"^"), < +00}. 

By [HI Proposition 2.2] the space G'5'BT/p(^, M™) coincide with {GSBVp{A,M.))"^, that is u 
. . . , u„,) G GS'BT/P(^, M™) if and only if u, G GSBVp{A, R) for every i = 1, . . . , m. 
The following compactness and lower semicontinuity result will be used in the following sections. 
For a proof, we refer to [2]. 

Theorem 2.1. Let A be an open and bounded subset o/R". Let g{x,u) : A x — > [0, 00] be a 

Borel function, lower semicontinuous in u and satisfying the condition 

lim g{x,u) = +00 for a.e. x € A. 

\u\ — >oo 

Let {uk)keN be a sequence in GSBVP{A;W^) such that 

limsup / \Vukix)\P dx + {S{uk)) + / gix,Uk{x)) dx < +00. 

k J A J A 

Then there exists a subsequence {uk^)heN o,nd a function u G GSBVP{A]W^) such that 

(2.1) Uk,^ u in measure, 

Vufe^ ^ Vu weakly in LP{A; M"^"). 

Moreover we have that 

n"--^ {S{u)) < liminf (S{ukJ) . 

h 

Let g g]1, +00 [ and let us set 

(2.2) GSBVP{A-W) GSBVP{A;R"') n L«(^;M'"). 
We say that Uk ^ u weakly in GSBVP{A;R"') if 

Uk ^ u in measure 

(2.3) Vufc ^ Vm weakly in LP{A; Af™^") 

Uk-^u weakly in L«(A;R'"). 

We will often use the following fact: li Uk ^ u weakly in GSBVP{A;W^) and F C A is such that 
Ti^~^{T) < +00 and S{uk) C F up to a set of 7i^~^-measure zero for all fc, then S{u) C F up to 
a set of 



F-convergence. Let us recall the definition of De Giorgi's T-convergence in metric spaces: we 
refer the reader to 11 for an exhaustive treatment of this subject. Let {X, d) be a metric space. 
We say that a sequence Fh X [—00, +00] F-converges to F : X [~oo, +00] (as h +00) if 
for all M G X we have 

(i) {T-liminf inequality) for every sequence (uh)/ieN converging to u in X, 

\im mi Fh{uh) >F{u); 

(ii) {T-limsup inequality) there exists a sequence {uh)heN converging to u in X, such that 

limsup Fh{uh) < F{u). 
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The function F is called the F-limit of (Ffi) (with respect to d), and we write F = F— lim/i F/j. 

We say that a family of functionals {F^} F-converges to as £ ^ if for every sequence Eh 
a.s h ^ +0O we have F— lini/j = F. 

The peculiarity of this type of convergence is its variational character explained in the following 
proposition. 

Proposition 2.2. Assume that the sequence {Fh)hefi T-converges to F and that there exists a 
compact set K C_ X such that for all h G N 

inf Fhiu) = inf Fhiu). 
ueK uex 

Then F admits a minimum on X , mix Fh —> minx F, and any limit point of any sequence {uh)heN 
such that 

lim (Fh{uh) - inf Fh{u)) = 

is a minimizer of F. 



Hausdorff metric on compact sets. Let A C be open and bounded, and let IC{A) be the 
set of all compact subsets of A. IC{A) can be endowed by the Hausdorff metric dn defined by 

duiKi, K2) max< sup dist(a;, i4r2), sup dist(j/, ifi) > , 
Ueifi yeK2 } 

with the conventions dist(x, 0) = diam(A) and sup0 = 0, so that dH{9,K) = if if = and 
dni^jK) = diam(A) if if 7^ 0. It turns out that IC{A) endowed with the Hausdorff metric is a 
compact space (see e.g. EI)- 

3. The quasistatic crack growth of Dal Maso-Francfort-Toader 

In this section we describe the quasistatic evolution of brittle fractures proposed in . They 
consider the case of n-dimensional finite elasticity, for an arbitrary n > 1, with a quasiconvex bulk 
energy and with prescribed boundary deformations and applied loads, depending on time. Since 
we are going to approximate the case n = 2, we prefer to introduce the model in this particular 
case. For more details, we refer the reader to jl3| . 

Let 17 be a bounded open set of with Lipschitz boundary and let fis be an open subset of 
n. Let Q d^l be closed in the relative topology, and let ■= \ On^- Let ds^ Q 
be closed in the relative topology and such that Qb ds^ = 9- In the model proposed in [T^ . 
^Ib represents the brittle part of Jl, and the part of the boundary on which the deformation 
is prescribed. Moreover the elastic body fl is supposed to be subject to surface forces acting on 

dsn. 

Admissible cracks and deformations. The set of admissible cracks is given by 

n(TiB;dNn) {F : F is rectifiable , F C (flis \ ^tv^?), n\T) < +00}. 

Here Ac B means that A C B up to a. set of Ti^-measure zero, and F rectifiable means that there 
exists a sequence (Mi) of C^-manifolds such that Fc Ui-^i- 1^ ^ is rectifiable, we can define 
normal vector fields to F in the following way: if F = IJ^ F^ with F^ C Ali and F^ n Fj = for 
i ^ j, given x G F^, we take vix) — VMi{x), where VMi{x) is a normal vector to the C^-manifold 
Mi at X. It turns out that two normal vector fields associated to different decompositions IJ^ F^ 
of F coincide up to the sign almost everywhere. 

Given a crack F, an admissible deformation is given by any function u G GSBV{Vl\ M'^) such 
that S{u) C F. 

The surface energy. The surface energy of a crack F is given by 
(3.1) j k{x,v{x))dV}{x), 
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where v is a. unit normal vector field on T. Here /c : x R'^ — > M is continuous, k{x, ■) is a norm 
in for all x e Qb and for all x <E Qb and v gM.'^ 

(3.2) K^Wl <k{x,:y) <K2W\, 

where Ki,K2 > 0. Notice that since k is even in the second variable, we have that the integral 
H3.1|l is independent of the orientation given to F, that is independent of the particular choice of 
the unit normal vector field v. 

The bulk energy. Let p > 1 be fixed. Given a deformation u e GSBVP{il;M.'^) the associated 
bulk energy is given by 

(3.3) >V(Vu) := / W{x,\/u{x))dx, 

Jn 

where W : fl x Af [0, +cxd) is a Caratheodory function satisfying 

(3.4) for every x G ft : W{x, ■) is quasiconvex and on Af^^^, 

(3.5) for every {x, e x M^""^ : a^\^\P - (x) < W{x, < afl^r + b^i^)- 

Here ag^,a]^ > 0, and bg^jb^ & L^{n) are nonnegative functions. Quasiconvexity of W means 
that for aU ^ G M^^^ and for all € C^ifl; M^) 

W{^) < I W{£. + Vip) dx. 

If we denote hy d^W : ^ x M^^^ ]VP^^ the partial derivative of W with respect to ^, it turns 
out that there exists a positive constant > and a nonnegative function S (fi), with 
p' ■= p/[p - 1), such that for all (x, G x Af^xs 

(3.6) \d^W{x,0\<o^W-'+h'^{^)- 

By 1(331) and (ESJ the functional W, defined for aU $ e LP(f]; Af2x2) by 

>V($) := / W[x,^{x))dx, 
Jn 

is of class Ci on LP(f7; M^""^), and its differential : LP{n; Af^xZ) ^ ^^'(f]. m^""^) is given by 

(a>V($),^')= / d^Wix,^ix))^{x)dx, e Lf(f^;Af2><2)^ 

Jo 

where (•, •) denotes the duality pairing between the spaces LP (f7;Af2x2) and LP{n; IVP""^). By 
and lEiSIl, there exist six positive constants > 0, a];^ > 0, > 0, > 0, > 0, 
>0 such that for every ^' e LP(f7; A/2x2) 

(3.7) \{dwi^),n < («rii'i>iir'+/52^)ii*iip- 

The body forces. Let q > 1 he fixed. The density of applied body forces per unit volume in the 
reference configuration relative to the deformation u at time t G [0,T] is given by dzF{t, x,u{x)). 
Here F : [0, T] x x ^ M is such that: 

for every z G : (<, x) — + F{t, x, z) is x £^ measurable on [0, T] x fi, 

for every (t, x) G [0, T] x 17 : z ^ F(t, x, z) belongs to C^(R2), 

and satisfies the following growth conditions 

(3.8) ao|z|«-6o(t,a;) < -F{t,x,z) < af + 6f (t, x), 

\d.F{t,x,z)\<a^\z\^-'+b^{t,x) 

for every (t,a::, z) G [0,r] x x E^, with > 0, af > and af" > 0, and where bQ,bf e 
C°([0,T];Li(r2)), 6f e CO([0,T];L9'(r2)) are nonnegative functions, with q' := q/{q - 1). 
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In order to deal with time variations, we assume also that for every {t, z) e [0, T] x 

F{t,x,z) = F{0,x,z) + I F{s,x,z)ds for a.e. a; G 
Jo 

dzF{t, X, z) = dzF{0, x,z) + / dzF{s, x, z) ds for a.e. x gQ, 

Jo 

where F : [0, T] x x ^ M is such that 

for aU z e : (t, x) ^ F{t, x, z) is x measurable on [0, T] x Q, 

for all (i, a;) G [0, T] X : ^ ^ F{t, x, z) is of class on , 

and satisfies the growth conditions 

\F{t,x,z)\<al{t)\z\^ + hl{t,x), 

\dzF{t, x,z)\< af {t)\z\'^-^ + 6f (i, x) 

for all (i,x,z) e [0,r] x O x R'^. Here 1 < g < 9, and of ,af e L^{[0,T]), 6f e ii([0,r];Li(O)), 
64 S L^([0,T];L^ ^Q)) are nonnegative functions with q' := 

Under the previous assumptions, for every t € [0, T] the functional 

(3.9) J='{t){u):= / F{t,x,u{x))dx, j^{t){u) := I F{t,x,u{x)) dx 

Jq Jn 

are well defined on L'(0;M^) and i'(0;M^) respectively. Moreover we have that J^{t) is of class 
on L«(0;M2)^ ^ith differential dJ^{t) : Li{Cl;R'^) L«'(0;M2) defined by 

{dJ^{t){u),v) = / dzF{t,x,u{x))v{x)dx, u,w e L''(0;M^), 

where (•, •) denotes now the duality pairing between (r2;IR^) and L«(f2;R^). !F{t) is on 
L«(0;M2) with differential defined by 

{dT{t){u),v) = / 92i^(t,a;,u(a;))w(x)dx, it, ?; G R^), 

where (•, •) denotes the duality pairing between L'i' {fl;M.'^) and L'i{fl;M.'^). For all u,v € L9(fi;IR2) 
and for all t G [0, T] we have 

jr(t)(u) = .F(0)(w) + / j'{s){u) ds, 
Jo 

(3.10) {dJ^{t){u),v) = {dJ^{0){u),v)+ f {dJ^{s){u),v)ds. 

Jo 

Moreover we have that for every t G [0, T] and for every u, v £ L^(f2; R") 

ll^ll^ -Po < -miu) < af\\u\\l + /3f , 

(3.11) \{dJ^{t){u),v)\ < {a^\\u\\l-'+(3^)\\v\\„ 

(3.12) \m{u)\<a^mu\\i + Pi{t), 

(3.13) I < K{t)Mr'+0fmv\u, 

where aQ > 0, af > 0, af^ > 0, I3q > 0, /3f^ > 0, /J^^ > are positive constants, and 

, af , , f3^ G L^([0,T]) are nonnegative fimctions. 

The surface forces. The density of the surface forces on Sgfi at time t under the deformation 
u is given by dzG{t, x, u{x)), where G : [0, T] x 9sf2 x ^ R is such that 

for every z gM."^ : {t, x) — > G{t, x, z) is x W^-measurable, 

for every {t,x) G [0,T] x 9sf2 : z G{t,x,z) belongs to C^(M2), 
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and satisfies the growth conditions 

- a^{t, x)\z\ - b^{t, x) < -G{t, X, z) < a'{\z\'' + bf{t, x), 
\dzGit,x,z)\ <a^\zr' + b^{t,x), 

for every {t,x,z) G [0,T] x 9^0 x M?. Here r is an exponent related to the trace operators on 

Sobolcv spaces: if p < 2, then we suppose that p < r < tt—, while if p > 2, we suppose p < r. 

Moreover af^ > 0, > are two nonnegative constants, and S L^{[0, T]; U (ds^)), b^, bf € 
C°{[0,T];L\ds^)), and b§ e C°{[0,T]; U' {ds^)) are nonnegative functions with r' := r/{r - 1) 
We assume that for every {t, z) G [0, T] x 

G{t,x, z) = G{0,x, z) + / G{s,x,z)ds for Ti^-&.e. x G ds^, 
Jo 

dzG{t,x, z) = dzG{0,x, z) + / dzG{s,x,z)ds for W^-a.e. a; G 9^0, 

JQ 

where G : [0,T] x ds^ x ^ K is such that 

for all z e : (t, x) G{t, x, z) is x H^-mcasurablc, 

for all {t,x) e [0,r] X ds^ : z G{t,x,z) belongs to C\R^), 

and satisfies the the growth conditions 

\Git,x,z)\<a^{tM- + b^{t,x), 

\dzGit, x,z)\< a2{t)\z\'-^ + b^{t, x) 

for all {t,x,z) G [0,r] x ds^ x R^. Here a§, af G L\[0,T]), bf G L\[0,T]; L\ds^)) and 
bf G L^([0,T];iy'' (ds^)) are nonnegative functions. 

By the previous assumptions, the following functionals on i''(9sf2;M^) 

(3.14) g{t){u):=( G{t,x,u{x))dn\x), g{t){u) := ( G{t,x,u{x))dH^{x) 

JdsO JdsO 

arc well defined. For every t G [0,T] we have that Q{t) is of class on L'^ {dsfl-jM."^) and its 
differential is given by 



{dg{t){u), v)= [ dzG{t, X, u{x))v{x) dH^ (x), u,v G L^dsfl; 



Uz^\L, x, uyxjjv\xj un. \x), u, f t jt"^^ 

where (■, •) denotes now the duality pairing between L'^ (9sfi;M^) and ^'"(Qsfi; R^). Moreover, 
g{t) is of class on L'^{dsi^; R^), and its differential is given by 



{dg{t)iu),v) = / dzG{t,x,uix))v{x)dn^{x) 

for all u,v G L^{ds^;M.'^)- Finally we have 

Q{t){u) = g{0){u) + I g{s){u)ds, 
Jo 

{dg{t){u),v) = {dg{o){u),v) + [\dg{s){u),v)ds, 

Jo 

for every u,v G L^{ds^;S?). 

Let Cts C i7 \ ilg be open with Lipschitz boundary, and such that ds^ C dfls', the trace 
operator from W^'P{^Isj^^) into L''(5f25;R^) is then compact, and so there exists a constant 
7s > such that 

(3.15) \\u\\r,ds(i < 7s(l|Vu||p,ns + l|u||j),ns) 
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for every u g W^^'^ifls', M^). By the previous assumptions, we have that there exist six nonnegative 
constants a^, af, af, , /3f, /3f and four nonnegative functions af, af, /Jg, G L^{[0,T]), 
such that 

-a^\\u\\r,9sn-P^ < -gmu) < a?h||;,asO 

(3.16) \{dg{t)iu),v)\ < {al\\u\\l;gl^+(3^)\\v\\r,asn, 

(3.17) \m{u)\<aimu\\:^,^^+piit), 

\{dg{t){u),v)\ < {ai{t)\\u\\l;gl^+p!m\v\kosn 
for every t e [0, T] and u,v e L''{ds^; M^). 

Configurations with finite energy. The deformations on the boundary do^ are given by (the 
traces of) functions g e W^^'P(17;R^) n L'?(f2;M^), where p,q are the exponents in H3.5(l and (|3.8|l 
respectively. Given a crack T G 7?,(riB; c^^rri) and a boundary deformation the set of admissible 
deformations with finite energy relative to (5, F) is defined by 

AD{g, V) {u 6 GSBVP{n] R^) : S{u) c T, u = g H^-a.e. on do^ \ T}, 

where we recall that 

and the equality u = g on dn^X^ is intended in the sense of traces (see Section 2]). 

Note that if u e G S EVP {n;R'^), then W{u) < +00 and \T{t){u)\ < +00 for all t € [0,T]. 
Moreover since T e 7^(^i3; ^Arf]) and S'(u) C T, we have that u e W^'PinsiR"^) D ^'(f^s;^^) so 
that Q{t){u) is well defined and \Q{t){u)\ < +cx) for all t G [0,T]. Notice that there exists always 
a deformation without crack which satisfies the boundary condition, namely the function g itself. 

The total energy. For every t G [0,r], the total energy relative to the configuration (u,r) with 
u G AD{g, r) is given by 

(3.18) £{t)iu,r) ■.= £''^{t){u)+£%T), 
where 

(3.19) S^'mu) ■■= W{u) - T{t){u) - g{t){u), 

and yy, T{t), g{t) and are defined in (|^ . ifT^ . and l|TT|l respectively. It turns out 

that there exist four constants > 0, af > 0, /3o > 0, /3f > such that 

(3.20) £^\t){u)>al,{\\Vu\\l+\\u\\l)~l3l 

£'\t){u) < af {\\Wu\\P + \\ur^ + + 
for every t G [0,T] and u G G5ByP(r2; K^). 

The time dependent boundary deformations. We will consider boundary deformations g{t) 
such that 

t g{t) G AG{[0,T];W^'P{n;R^) n L'^{n;R'^)), 

so that 

t g{t) G L^([0,T];W^i'f(f7;R2) ni«(f7;R2)), 

and 

V5(t) Gii([0,T];LP(r!;Af2^2))_ 

The existence result. Let Fq G TZ{flB,dN^) be a preexisting crack. The next Theorem proved 
in establishes the existence of a quasistatic evolution with preexisting crack To. 

Theorem 3.1. Let Tq G TZ{Qb] disiQ) be a preexisting crack. Then there exists a quasistatic 
evolution with preexisting crack Tq and boundary deformation g{t), i.e., there exists a function 
t {u{t),r{t)) from [0,T] to GSBVP{n;R'^) x TZ{nB;dN^) with the following properties: 
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(a) (u(0),r(0)) is such that 

£(0)(m(0), r(0)) = min{f (0)(z;, T) : t- e AD[g{Q),T), Tq C T}; 

(b) u{t) G AD{g{t),V{t)) for all t £ [0,T]; 

(c) irreversibility: Tq ClT{s)(ZT{t) whenever < s < t < T ; 

(d) static equilibrium: for all t G [0, T] 

£{t){u{t),T{t)) = mm{£{t){v, T) : v e AD{g{t),T), T{t) C T}; 

(e) nondissipativity: the function t E(t) :~ £(t)(u(t),T{t)) is absolutely continuous on 
[0,T], and for a.e. t G [0,T] 

(3.21) m = (dW{Vuit)),Vg{t)) - {dJ'{t){u{t)),git)) " 

-{dg{t)(u(t)),g{t))-Q{t){u{t)). 

The next theorem gives a compactness and lower semicontinuity result with respect to weak 
convergence in GSBVP{fl,M.'^) which will be often used in the next sections. 

Theorem 3.2. Let tk G [0,T] with tk t, and let (uk) C GSBVP{n;R^), C g]0;+oo[ such that 
S(uk) C flB and 

r\tk){uk)+£'iSiuk)) < C, 

where £^^ and £^ are defined as in 13.19() and . Then there exists a subsequence (ufc^)heN 
converging to some u weakly in GSBV^{Q,]^) such that S{u) C VLb, 

r\t){u)<\ijnmi£^\tkj{uk,). 

and 

£%S{u)) < \immi£%S{ukJ). 

h^oo 

Proof. By (|3.20|1 and (|3.2|1 . we have that there exists C g]0, +oo[ such that 

\\vuk\\; + \\uk\\i + n\s{uk))<c'. 

Then we can apply Theorem 12.11 with g(x,Uk) = \'Uk\'^ , obtaining a subsequence (ufc,J/igN and 
u G GSBV'''{^\R^) such that (|2.1() holds: in particular we may assume that — > u pointwise 
a.e.. We have u^^^ u strongly in L^{Vl\M?), and by Fatou's Lemma we have that u G L'(ri;R^) 
so that u G GSBVP(n\m.^). We conclude u weakly in G5By/(r2; R^). By 2, Theorem 

3.7] we have that 

<liminf£^(5KJ), 

h 

by we have that 

/ VK(x, VM)da; < liminf / W{x,Vuk^)dx, 
Jn ^ Jn 

and by Fatou's Lemma (in the limsup version) we have 

limsup / F{tk^,x,Uk^{x))dx < / F{t,x,u{x)) dx. 
h Jn Jn 

Since (w,,J|f^^ is bounded in W^^P{ns;^'^)r\Li{fls;M.'^), and the trace operator from W^i'P(f]s; K^) 
into {ils;M.'^) is compact, we get 

iimgitk,)iukj^g{t)iu), 

h 

and so the proof is thus concluded. □ 
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Figure 1. 



4. The finite element space and the transfer of jump 

Let i7 C R2 be a polygonal set and let us fix two positive constants < Ci < C2 < +00. By 
a regular triangulation of of size e we intend a finite family of (closed) triangles Tj such that 
fl = [j-Ti, Ti Ci Tj is either empty or equal to a common edge or to a common vertex, and each 
Tj contains a ball of diameter ciS and is contained in a ball of diameter C2£. 

We indicate by 7ls{^) the family of all regular triangulations of fl of size e. It turns out that 
there exist < i?! < ■(92 < tt such that for all triangle T belonging to a regular triangulation 
T e TZs{Q), the inner angles of T are between t?i and '&2- Moreover, every edge of T has length 
greater than Ci£ and lower than C2£. 

Let e > 0, Rg e 7?.e(0), and let a g]0, i[. Let us consider a triangulation T nested in 
obtained dividing each triangle T e Re into four triangles taking over every edge [a;, y] of T a knot 
z which satisfies 

z = tx + {1 — t)y, t € [a,l — a]. 

We will call these vertices adaptive, the triangles obtained gluing these points adaptive triangles, 
and their edges adaptive edges. We denote by 7^^a(r2) the set of triangulations T constructed in 
this way. Note that for all T G 7^_o(f2) there exist < < C2 < +00 such that every Tj e T 
contains a ball of diameter c"e and is contained in a ball of diameter CjE. Then there exist 
< < ^2 ^ such that for all T belonging to T £ %^a{^), the inner angles of T are between 
•&1 and '&2- Moreover, every edge of T has length greater than c°£ and lower than C2£. 

Prom now on for all £ > we fix R^ € TZe{fl). We suppose that the brittle part Qb and 
the region Qs introduced before for the model of quasistatic growth of fractures arc composed of 
triangles of Rg for all s. Moreover we suppose that and ds^ are composed of edges of Rg 
for all £ up to a finite number of points. Finally, in order to deal with the deformation at the 
boundary, it will be useful to consider iljj polygonal such that CId D f2 = 0, and dCtn dCl = 
up to a finite number of points. We set 

(4.1) n' ■.= nunDLidD^, 

and we suppose that the regular triangulation Rg can be extended to a regular triangulation of 
so that every triangulation T in 7^^o(f2) can be extended to a triangulation of T^^ai^') considering 
the middle points of the new edges as adaptive vertices: we indicate these extended triangulation 
with the same symbol T. 

We consider the following discontinuous finite element space. Wc indicate by Ae,ai^',^'^) the 
set of all M : i7 ^ such that there exists a triangulation T(?i) £ 7^^a(i7) nested in R^ with u 
affine on every triangle T £ T(m). For every u £ Ae,a{^',^'^), wc indicate by S{u) the family of 
edges of T(u) inside across which u is discontinuous. Notice that u £ SBV{^; R^) and that the 
notation is consistent with the usual one employed in the theory of functions of bounded variation. 
Let us set 



(4.2) A7-£(ri;K^) :={-«: O ^ : u is continuous and affine on each triangle T £ Re}. 
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The projection TTaiP) The interpolating curve S^.a 




Figure 2. 

The discretization of the problem will be carried out using the space 

(4.3) Afjrt; M^) {u g A^A^; : Siu) C Ub}. 
Given any g e AF£(f7;R2), ^^^^^ ^ ^ Af^^^R"^) let 

(4.4) Sl,{u) := {x € don : u{x) + g{x)}, 

that is S\)(u) denotes the set of edges of at which the boundary condition is not satisfied. 
For every u £ ^^£,(17; M^), let us also set 

(4.5) S^{u):^S(u)\JS%(u). 

An essential tool in the approximation result of this paper is Proposition 14 . 21 which generalizes 
the piecewise afHne transfer of jump |17[ Proposition 5.1] to the case of vector valued functions 
with bulk energy E'^^ and surface energy f of the form (|3.19l) and H3.1II respectively. 

In order to deal with the surface energy f we will need the following geometric construction. 
Let 5' C be a segment and let us suppose that S intersects the edges of Re at most in one point 
for all £ > 0. Let a g]0, \ \, and let P = S' fl C, where C, = \x,y\ is an edge of R^: we indicate 
with TTa{P) the projection of P on the segment {tx + (1 — t)y : t E [a, 1 — a]}. The interpolating 
curve Se,a of S in R^ with parameter a is given connecting all the 7ra(P)'s belonging to the same 
triangle of R^ (see Figure 2). 



Lemma 4.1. Under the previous assumptions, there exists a function r]{a) independent of S with 
r]{a) as a ^ such that 



limsup \£'{Se.a) - £"{S)\ < r){a)£'{S), 

e->0 



where is defined in (|3.1|l . 



Proof. By ()3.2|l . we have that there exist w and > such that for all xi^X2 G ^ and |i^i| = 

k2| = 1 ' 

\k{xi,vi) - k{x2, V2)\ < (^i\xi - X2I) + K^li^i - 1^2!, 
where ui : ]0, +oo[^]0, +oo[ is a decreasing function such that io{s) ^ as s — > 0. Let T e R^ be 
such that T n 5 ^ 0, and let us choose xt T D S and xfp'^ e T n 5^ q. Let C2 > denote the 
characteristic constant of R^ such that every T G R^ is contained in a ball of diameter C2S. Then 
we have 



/ k{x, ly^"-) dn^ - / fc(x, j/t) dH^ 

/ k{x'^r^°' ^vfp°') dTi} — I k{xT T vt) dTi^ 
Jsr„nT JsnT 



+ uj{c2e)n\Se,a n T) + uj{c2e)n\s n t) 



< \k{x'T"',VT'')n\Se,anT) - k{xT,i'T)'H^iS nT)\ + Uj{c2e) [H^ {Se,a ^ T) + H^S H T)] , 
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where t'y", i^t are the (constant) normal to S^^a H T and S HT respectively. We have 

|fc(a:^",4'")H^(5e,a nT) - k{xT,iyT)n^{snT)\ 

< fc(a;^'',4'°) \H\Se,ariT)~H^{SnT)\ + |/c(x^°, 4''') - k{xT,iyT)\n\S DT) 

< K2 |7^^(S'e,anT) -H^(S'nT)| +tj(|a;^" -a;T|)H^(5nr) +ii:3|4'" -z^T|H^(5nT), 

where K2 is defined in %\.2\ . We are now ready to conclude: in fact, following |19l Lemma 5.2.2], 
we can choose the orientation of ° in such a way that 

14'" - z^tIH^STiT) <D2ae, | H^S'e.a n T) - (5 n r)| <Diae, 

with Di,D2 > independent of T, e,a. Then, summing up the preceding inequalities, recalling 
that the number of triangles of Rg intersecting S is less than ce~^TC^{S) for e small enough, with 
c independent of S and e (see for example ^| Lemma 2.5]), we obtain 

limsup \£%Se.a) - £'{S)\ < p{a)H^{S), 

where p{a) :— c{K2Di + K3D2)a. In view of H3.2|l . we conclude that 

limsup \£''{Se,a) - £%S)\ < K^^p{a)£%S), 

and so the proof is concluded choosing 77(a) := K^^p{a). □ 

For aU u £ GSBVP{n;M.'^) and for aU g £ W^'P{n;M.'^) n L«(17;M2)^ let us set 
(4.6) S^iu) S{u) U {x £ do^ : u{x) ^ 5(2;)}, 

where the inequality is intended in the sense of traces. We are now in a position to state the 
piecewise affine transfer of jump proposition in our setting. 

Proposition 4.2. Let a e]0, i[, and for all i — 1, . . . ,ni let 

ul£Afjn;R^), u' £GSBVP{n;R^) 

be such that 

ul u' weakly m GSBVP[n- M^). 
Let moreover gl,he £ ATe{^;R'^) , 9\h£ W'^^P{n;M?) n L''{n;M?) he such that 
gl g\ h^^h strongly m W^'P{n; M^) n L'(17; M^). 
Then for every v £ G'5ByP(f7; R^) with S{v)cQb, there exists £ ^fa(f^;M^) such that 

Vwe -> S7v strongly m LP{VL\ M^'''^), 
^ V strongly in L'^(ri;M^), 

and such that 

limsupf^ (^S'^^iv,)\Qs^Hul)^ <p{a)£' (^S''iv)\[j^S<'\u^)^ , 

where p,{a) depends only on a, /i(a) -^1 as a — > 0, and £^ is defined in In particular for 

all t £ [0, T] and for all t^ t we have 

£^\t,){v,) ^ £^\t){v), 

where £'^^ is defined in (|3.19|) . 

The proof of Proposition 14.21 can be obtained from that of [IT, Proposition 5.1] taking into 
account the following modifications. We can consider v scalar valued since vector valued maps can 
be easily dealt componentwise. Even if the surface energy is of the form Ij^.l!) , we can still restrict 
ourselves to the case in which v has piecewise linear jumps outside a suitable neighborhood of 
yJiLi i"^^) using the density result of 10 . In order to approximate the piecewise linear jumps, 
we use Lemma |4.1I Finally the fact that p ^ 2 prevents us from considering the piecewise jumps 
as union of disjoint segments: we overcome this difficulty choosing = in the regular triangles 
which contain the intersection points, and then interpolating v outside as in (171 Proposition 5.1]. 
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5. Preexisting cracks and their approximation 

In Sectional wc will need to approximate the surface energy of a given preexisting crack F*^. 
We take the initial crack in the class 

(5.1) T{n) {r cUb : H\r) < +00, T = S''{z) 

for some h e W^''''{n;R^) D L''{n;M.^) and z e G5SV'P(f7; M^)}. 

Notice that it is not restrictive to assume h = 0. We take as discretization of r(ri) the following 
class 

(5.2) Te,a{^) := {rcHs : ^^^(r) < +oo, T = 5°(z) for some z G Af^jn;R^)}. 
We have the following approximation result. 

Proposition 5.1. Let r° e Then for every e > and a e]0, i[ there exists r° ^ e F^^ai^) 

such that 

lim r(r°,j = r(r°), 

e,a — >{) 

where S'^ is defined in 

Moreover let G AFe{^;'^^), 9 e H/i^P(r2; M^) n L9(f7; R^) be such that as s ^ 

Qe^g strongly m R^) n M^), 

and let us consider 

_ /f'^'(0)(«)+f«(59'(«)\rOJ z/«G^^Jf};R2), 

foo otherwise in L^{Vl\M?), 



^^^^ ^ |£^'(0)(«)+f^(59(«)\r0) ifveGSBVP{n;R^),S{v)cnB, 
1 +CXD otherwise in L^{fl;M.'^), 

where S^^ is defined in (|3.19(l . Then the family {F^A ^-converges to F in the strong topology of 
{n;M.'^) as e and a ^ 0. 

Proof Let us consider T" £ T{n) with r° = 5°(z) for some z £ GSBVPifl;^?). Then by 
Proposition 14 . 21 for every e > and a G (0, i), there exists z^^a G -4e.a(f^; IR^) such that for e — * 
and for all a 

VS^^a Vz Strongly in LP{n; M^""^), 
Ze.a z strongly in L^(r2; 



and 



limsupf (5''(5e,a)) < iJL{a)£'{S''{z)) 

e-*0 

with ^{a) ^ 1 as a ^ 0, where £^ is defined in (|3.1|l . Let \ 0, and let \ be such that for 
all e < Ei 

£'{S"{~z,^a^)) < pi{a,)£'{S"{z)) + a,, 

and 
Setting 



we have that 



and 



lim Vze.a = Vz strongly in i^(fi; Af^' 
lim a = -2 strongly in M^), 



limsupf«(5"(z,,,)) <r(5"(z)). 

e,a—^0 
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Since by Theorem 13. 21 we have ^^(^"(zj q)) < hminf^ a^o ^-*('5'°(^e,a)), we conclude that 

hm 5^(^0(z,,,))=£^(5°(z)). 

£,a— »0 

Let us set for every £, a 
We have that 

um r(r«,j = f^(r°). 

Let us come to the second part of the proof. Let us consider (£„, a„)„gN such that £„ — s- and 
a„ — > 0. If we prove that (Fe^^a„)neN F-converges to F in the strong topology of L^(ri;M^), the 
proposition is proved since the sequence is arbitrary. Since we can reason up to subsequences, it 
is not restrictive to assume a„ \ 0. 

Let us start with the F-limsup inequality considering v e GSBV^{Q,]A^), with S{v) C {Ig. For 
any n fixed, by Proposition 14. 21 there exists Ve^a„ G -^fa„(^!^'^) such that for e ^ 

Vue,Q„ Vv strongly in LP{Vl] M'^'^'^), 

V£,a„ V strongly in L'^(fi;M^), 

and such that 

limsupr(5»-(«e.aJ\F2, )<^iia,,)£'iS3^iv)\T") 
with /i(a) 1 as a ^ 0. For every m G N let e™ be such that for all e < e™ 
£'{S^^{i>e,aJ \ ri^J < fi{a^)£'{SS{v) \ F") + a™, 

and 

l|VCe^<j,„ - V^;||ip(0;M2x2) < a^, l|we,a„ - v\\Li{a-R^) < am- 
We can assume e™ \ 0. Setting 



_m+l 



we have that 



and 



= Vv strongly in LP{n; M^""^), 
limWe^^Q^ — V strongly in L"^ {Q;'R'^), 



limsupr(5«'(i;e„,aj \F^_J < £%SSiv)\T°). 

n 

Then we get 

limsupF,„,,„(^;,„,aJ < limsupr'(0)K„,a„) + limsup^ (^f- K„,a„) \ F^^^.J 

n n n 

<£'^'iO)iv)+£%S3{v)\r°)^Fiv), 

so that the F-limsup inequality holds. 

Let us come to the F-liminf inequality. Let £ L^(il;R^) be such that w„ v strongly 

in Li(f];M2) and liminf„ Fe„, < +oo. By Theorem H!^ we have v e GSBVP{n;R^) with 
S{v) cTIb and 

£''\0)iv) < liminf£:"'(0)(w„). 

n 

Let us consider fl' defined in Let us extend g^^ and g to W^'P{n';R'^)nLi{n';R'^) in such a 

way that g strongly in W^'P{n'; M^) n L'^(ri'; M^), and let us also extend v to fi' setting 

?;„ = g^^ and u = g on fi^). We indicate these extensions with w„ and w respectively. Let us also 
set 2;£„,a„ ^ z = on where ZE„,a„ and z are such that F^ ^ = z^^ a^ and F° = S'(z). We 
indicate these extension by C£„,a„ and C respectively. Then for every 77 > we have by Theorem 

EH 

£'iSiw + j?C)) < liminf r (5K + r,Ce„,aJ). 
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Since for a.e. 77 > we have S{'w + j]C) — S{w) U S{Q and S{wn + JyCen.o™) = S{wn) U 'S'(Ce„_o„), 
we deduce that 

n n, n 

Since by assumption £^{T^^ — > .^^^(r'^), we conclude that 

r(5f(«)\r") <hmmfr(5f-K)\r° a )■ 

We deduce that 

s^'mv) + («) \ ro) < Hminf [r'(o)K) + ^ (5^- («„) \ ,J] 

that is 

F{v) < hminf i^.^.a„(u„). 

n 

The F-hminf inequahty holds, and so the proof is concluded. □ 

6. The discontinuous finite element approximation 

In this section we construct a discrete approximation of the quasistatic evolution of brittle 
fractures proposed in ^31 a-nd described in the Preliminaries: the discretization is done both in 
space and time. Let us consider 

ge e W^'\[Q,T];W^'P{n;R.'^) n L''{n;m.'^)), g,{t) e ATe{^;^^) for aU t G [0,T], 

where ^JFg(0;]R^) is defined in (|4.2|l . Let (5 > 0, and let Ns be the largest integer such that 
d{Ns - 1) < T; we set tf := iS ior < i < Ns - 1, := T and gl'^ := ge{4). Let G T^^a{^) 
be a preexisting crack in il, where r£^Q(ri) is defined in 1)5. 2|l . 

Proposition 6.1. Let e > 0, a g]0, ^[ and 5 > Q he fixed. Then for all i — 0, ... , Ns there exists 
""e.i £ -4fa(^;I^^) t/iat, setting 



r=0 

iB ('o.^^>2^ 



we /laue for all v G -4^a(ri; M ) 

(6.1) ^-^'W^",) + 5^ \ r°) < £^'(0)(7;) + £^ [s^'-\v) \ F") , 
and for I < i < Ns 

(6.2) S''(tt)iv^^ + £^ (5^"' («f:l) \ Ff:r^) < ^-^'Itf )(^) + (^;) \ rf;r') ■ 
Proo/. Let be a minimum of the following problem 

(6.3) min {£^\Q){n) + {S^"\u)\T^)] . 

We set Ff'" := F" U " (wf'J^). Recursively, supposing to have constructed u^^'^~^ and Ff'^"-'^, let 
wf'^ be a minimum for 



(6.4) min {£^\t^){u) + £:^(5''e"(«) \ Ff--^)} 



We set Ff;j, Ff-j,"! U ^f'"' (wf;*^). It is clear by construction that (jHU and hold. 

Let us prove that problem H6.4() admits a solution, problem (|6.3(l being similar. Let (un)nGN be 
a minimizing sequence for problem H6.4|l : since g^'^ is an admissible test function, we deduce that 
for n large 

£'^'{tf){u,,)+r (^^^^■'(^„) \Ff;r^) < e^'it'Mn + I. 

By the lower estimate on the elastic energy (I3.20f) . we deduce that for n large 



(6.5) 4 {\\^uJP + \\u„r^)+r (^^^'■'(w„) \Ff:r^) < s^'itfMn + 1 + /3, 
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Let US indicate by T^,. . . the triangles of T(u„). Up to a subsequence, there exists T = 
{T^,...,T''} e %.a{^) such that for alH = 1, . . . , fc we have in the Hausdorff metric (see 

Section El for a precise definition). Let us consider £ T, and let T* be contained in the interior 
of T*. For n large enough, is contained in the interior of T^; moreover is affine and by 

()6.5|l we have |Vu„|'' dx + H^nH^ooff ,.^2) < C, with C independent of n. We deduce that there 
exists a function affine on T* such that up to a subsequence m„ — > it uniformly on T*. Since 
is arbitrary, it turns out that is actually defined on and 

' n 

where denotes the restriction of £'^''{tf) to the maps defined on T'. Let u E Ae,a{^',^'^) 

be such that u = on T* for every i — I, . . . , k: we have 

< liminfr'(if)K). 

n 

On the other hand, it is easy to see that S^^^'' (u) is contained in the Hausdorff limit of S'^=''(u„), 
and so It G ^^„(f2;R^); moreover we deduce 

£^ [ss''\u) \ rf:r^) < limmf (s^'^u^) \ rf:r') • 

We conclude that m is a minimum point for the problem 1)6. 4|l . so that the proof is concluded. □ 



The following estimate on the total energy is essential in order to study the asymptotic behavior 

iait) :=«f;Lfor all <f < t < tf+i 



of the discrete evolution as (5 — * 0, £ — > and a ^ 0. Let us set wf ■— wf'a for all < t < t^,-^ 



andz = 0,...,iV5-l, <jr)-7.f;^^ 
Proposition 6.2. For all < j < i < Ns we have 

(6.6) f(tf)(«f:l,rf:l) < £(t|)(«fi,rfi) + [ ' (dW{S74jr)),S7g,ir)}dr 

' Hr){uUr)) dr - [' (a.F(r)(«L(r)), ^.(r)) dr 



where ^ — > as 5 ^ Q uniformly in e and a. 

Proof. By the minimality property (|6.2|) . comparing u*' \ with uf'^-i - gf'^-i + 5^'* we get 

(6.7) >v(v«f;i) - .F(tf )«;) - ^(if )«;) + S\S^' («f;l) \ rf_i) 

-e(tf)(«f;r'-fff''-^+5f^). 

We have 

(6.8) W(Vwf;r' - V.g,*''-i + Vg'/) = W(V<r') 

+ (aw(v<r' + ^f;r'(vgf'^ - V5f^*-^)), v^f^^ - v^f^'-^) 

,i5,i-l\ I / /ai/M/v7„,(5 I „.i5 



>V(V<r')+ (a>V(V<,(r)+<,(r)),Vg,(r))dr, 



where iJff^-i g]0, 1[ and vI^{t) i^?;jr^(V5f'' " ^9e'^^) for aU r G 
Similarly we obtain 

(6.9) .F(if)«r' - + 9^') - + / ' (aj-(tf)«Jr) + <,(r)), 5e(r)) dr. 
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and 



(6.10) g{ttM':-'^gt''-'+9tn = g{ttM%-') + / {dg{ttM^jT)+zUr)),g,{r))dT, 



where w^Jt) Af;rHffe'' " ffs'^"'), 4jr) := z/f.^Hsf" " ffe'^"') for all r G [tU,t'A. and 

Af;r\<'r' e]0,i[. 

Since by ifTTUIl we have for r e if [ 



(9.F(s)«,(r)+<,(r)),.g,(r))ds 



we get by 1)3.1311 



(6.11) |(a.F(t,*)(<Jr)+<Jr)),<7e(r))-(aJ-(T)(<,(r)+<,(r)),g,(r))| 

< [ ' \{dHs){uUr)+wUT)),g,{T))\ds 



< 



5.€,a II 



where 



5, 6. a 

It ■ 



max I hf;r' + Af;r'(fff"-5f'*-')ll|-i /_ < (s) + / /34^(s)ds 



l<i<Afi 



Similarly we obtain 
(6.12) 

where 



From (EHIl, taking into account (jnHJl, (ES, (EHH), JHHSl, we have 



(6.13) £(tf)(uf;;,rf;;) < £(tf_i)«r\ rf;r') + / (aW(V«f,„(r) +<,(r)), V5e(T))dr 
j-(T)«,(r))dr- /' (a.F(T)«,(r)+<,(r)),g,(r))dT 
a(r)(ul,(T))dr- /" (9(7(T)(uL(r)+zLW),ffs(r))dT 

+ 7^^'" / ll5e(T)||gC;T + 7^'''''' / \\ge(T)\r,ds^dT. 
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Taking now < j < i < N^, summing in Ht).13|l form tj to tf, we obtain 
(6.14) rffj < rfi) + / ' {dWiVuUr) + < Jr)), V5e(r)) dr 

HtMAt)) dr - r {dT{T){4jT) + <Jr)),5,(r)) dr 



G{r)i4jT))dT^ / {dg{T)i4jT) + zUr)),gAr))dT 

+ 7^^'" / UeiT^qdr + jg"-" / \\gAT)\\r,dsndT. 



Setting 

(6.15) /' |(9W(V<,(t) + < Jr)), V<7e(r)) - (aW(V<,(r)), Vge(r))| dr 

Jo 

+ f \{dnr){<.a{r)+wUT)),Ur)) (5.F(r)« Jr)), 5.(r))| dr 
Jo 

+ /' |(ag(r)(^.f,„(r) +zlJr)),5e(r)) - (ag(r)(^.f,„(r)), g,(T))| dr 
Jo 

+ 1^/'" \\9e{T)\\qdT + -fg-^''' \\g,{T)\\r,dsndT, 
Jo JO 

from (|6.14() we formally obtain (|6.6|l . Let us prove that e^g^^, ^ ^ as 5 ^ uniformly in e and a. 
By (|6.2() . comparing uf;^ with gf'^ and taking into account (|3.20(l . we get for alH = 1, . . . , Ns, 

l|V<L||p+hffJ|,<C', 

where 

C':=^. max (f'='(tf )(g,^'') + /3o') • 
i=o,...,Afi 

Since is Lipschitz, there exists Ks > depending only on p, q such that 

||w||p.ns < Ksi\\Vu\\p^ns + \\u\\q,ns) 
for all u e W^'P{ns;M.'^) n L9(175;R2). Taking into account l|TTK|l . we obtain 

||<'all.,asO < C" 

for some C" independent of (5. Since g, e W^i'i([0, T]; M^i'P(r2; R^) n L'i{n;R^)), we obtain that 
for all r e [0, T] as (5 ^ 

^'L(^) ^ strongly in LP{n; M^""^), 

wt.ai^) strongly in i«(f7;R2), 

4.ai^) strongly in L''{dsn;R'^). 

Moreover 7^'^'" ^ and jg^'"' ^ as (5 ^ 0. Finally, by ^21 Lemma 4.9], as (5 ^ we have that 
for all r G [0, T] 

|(9>V(V^.f,,(r) +<Jr)), Vg,(T)) - (aW(V<,(r)), Vge(r))| ^ 0, 
|(5^(r)(«f,„(r) +<,(r)),g,(r)) - (a^(r)(^f,Jr)),5e(r))| ^ 0, 

l(5e(r)(^f,„(r) +z,^(r)),<7e(r)) - (9a(r)«Jr)),5e(r))| ^ 0, 

uniformly in e,a. By the Dominated Convergence Theorem, we conclude that ^ as (5 — > 
uniformly in e and a, and the proof is finished. □ 
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7. The approximation result 

In this section we study the asymptotic behavior of the discrete evolution obtained in Sectional 
Let us consider a given initial crack r° £ r(ri) where T{^1) is defined as in H5.1|l . and a boundary 
deformation g G W^^^ {[0,T],W^'P{n;R^) n L''{n;R^)) . Let G T^^ai^) be an approximation 
of F*^ in the sense of Proposition 15. 11 and let us consider 

5e e M/i'i([0,r],T^i'''(17;R2)ni«(17;R2)), 

such that 

ge{t) G AFe{n;R^) for all t G [0,T], 

and such that for e ^ 

g, g strongly in W^'\[0, T], W^^P{n; E^) n M^)). 

Let 

{«Lrf;;), * = o,...,iva 

be the discrete evolution relative to the initial crack P^ ^ and boundary data g^ given by Proposition 
16.11 We make the following piecewise constant interpolation in time: 

(7.1) ^4^,{t):=v^^, TiM:=ri% gl{t) -.^ g^) tor tf < t < tl„ 

z = 0, . . . ,iV5 - 1, and uijT) := uf;^, TijT) := Pf;^^ (T) g,{T). 
By Proposition iniS for aU v G yl^^(f};K2) we have 

f ^'(o)«jo)) + (5«=^(«)(^f,,(o)) \ p° J < £^'(o)H + £^ [s^^'^^Hv) \ r?, J , 

and for aU t G and for aU v G Af_^{n;R'^) 

(7.2) f^'(if)«Jt)) < £'='(^f)(^;) (^^^'^ («) \ Pf^Jt)) . 

Here f*^' and £^ are defined in H3.19|l and (|3.1|l respectively. Finally for all < s < t < 1 we have 

(7.3) £:(tf)«Ji),rlJ0)<^(5-)«a(s),rlJs))+ [ \dW{VuijT)),Vg,{T))dr 

J 

^(r)«a(T)) - [\dT{r){uijr)),g,{T))dr 



Q{r){uia{r)) ~ / (5g(r)«,(r)),<?,(r)}dr + ef,,, 

where s| < s < s^^^ and tf < t < if^^, ef ^ is defined as in (|6.15|) . and f(t)(M,F) is as in (|3.18|) . 
Recall that e* ^ — > as (5 ^ uniformly in e, a. 

Comparing wf^(t) with (i) by ^T^, and in view of ^1^, |TTT|l . |TT^ . (|XT7j) . (|^ 

and H3.2|l . by (|7.3|) with s = we deduce that there exists C" g]0, +oo[ such that for all t, 6, e and 
a 

(7.4) II v<Jt)||, + Ikl JOII, + w'(rL(0) < C. 

By the time dependence of Oi in view of (|7.4() . by (|7.2|l and H7.3|l we have that there exists 
o* — > as (5, e ^ uniformly in a such that for all t G [0, T] and for all v G ^(f^; K^) 

(7.5) f^'(i)«a(i)) < f^'(i)(^^) 
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and for all < s < t < T 

(7.6) £m<ait), Tijt)) <£{s)iul,{s), Tijs)) + j\dW{VuijT)),Vg,{T)) dr 

t nt 



J S J S 

f Q{r){ula{r))~ [\dg{r){4jr)),g,ir))dr + ol,. 



Inequality (|7.4() gives a natural precompactness of (uf ^{t)) in GSBV^{n;M.'^). The main result 
of the paper is the following. 

Theorem 7.1. Let 6 > 0, e > 0, a e]0, i[, and let {t ^ rf,a(t)) : t e [0,r]} be the 

discrete evolution given by 1)7.1(1 relative to the initial crack ^ and the boundary data g^. Then 
there exist a quasistatic evolution {t —^ {u{t),T{t))} in the sense of Theorem VS . 1\ and sequences 
Sn — > 0, £„ — > 0, a„ — » 0, such that setting u„{t) a„(^) '^'^'^ ^n{t) '■— T*^ a„(^)j /o'' 

t e [0, T] the following facts hold. 

(a) For every t G [0,T], (w„(t))„gN is weakly precompact in GSBVP{rt;'R.'^), and every accu- 
mulation point u(t) is such that S'^*^*)(ti(i)) CT{t), 

(7.7) S''{t){u{t)) < £''{t){v)+£' (5s(*)(«) \r(t)) 
for all V e GSBVP{n;M.'^) with S'(w) cUb, and 

£-\t){u^{t)) ^ £-\t){u{t)). 
Moreover there exists a subsequence of (dm £m o,n)neN (depending ont) such that 
Unit) ^ u(t) weakly m GSBV.f (^2; R^). 

(b) For every t £ [0, T] we have 

(7.8) £{t)iu„{t), r„(i)) ^ £{t){u{t), T{t)); 

more precisely elastic and surface energies converge separately, that is 

(7.9) £'^'{t){un{t)) -> £'"{t){u{t)), £^(r„(t)) -> £'im). 

For the proof of Theorem l7. II we need two preliminary steps. First of all, we fix a and study the 
asymptotic for 5, £ — >• (Lemma I7.2|l , and then we let a — > using a diagonal argument (Lemma 

rrti . 

Lemma 7.2. Let a be fixed, t € [0,r], and let (5„ and £„ —>■ 0. There exists Ta{t) G 
TZ{CIb', On^) o,nd a subsequence of {5n,£n)n&i (which we denote with the same symbol), such that 
the following facts hold: 

(a) ifwn e Af^^^a.{^;R'^) is such that S^'"^*^{wn) C Ff^ „(t) and 

Wn^w weakly in GSBVP{n; R^), 

then we have 

S3^'\w)cTa{t); 

(b) there exists fJ.{a) with fj.{a) I as a —> such that for every accumulation point Ua{t) of 
("f:,a(*))«6N for the weak convergence m GSBVP{n-m?) and for all v G GSBV^in-R^) 
with S{v) C f^B; we have 

(7.10) £''\t){ua{t)) < £'\t){v) + ^Ji{a)^ (S<^^'\v) \ Va{t)) ; 
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moreover 

(7.11) \\Tn£^\t){ull^M) = £'\t){ua{t)); 
(c) we have 

n " ' 

Proof. We now perform a variant of |131 Theorem 4.7]. Let {'Pk)k&\ C Li(r2;R^) be dense in 
L^(ri;R^). For every ipk and for every m £ N, let w^'^(i) be a minimum of the problem 

min{||Vw||p+ ||i;||q + m||i;-(^fe||i : w £ F"}, 

where 

Since by \7A\ we have 'H^(rf^ < C", by Theorem l2.1l there exists a subsequence of ((5„, £„)„gN 
(which we denote with the same symbol) such that vj,''^{t) weakly converges to some v'^^{t) e 
GSBVP{n; M^) as n ^ +00 for all k,meN. We set 

(7.12) r„(i):=|j5^'W«„(i)). 

k,m 

Let us see that Ta{t) satisfies all the properties of the lemma. Clearly Tait) e TZ{flB- On^) and 
point (c) is a consequence of Theorem 12. II In particular by H7.4|l we have that 

(7.13) n\Ta{t))<C'. 

Let us come to point (a). Let w„ e ^f^_„(il;M2) be such that 5"^=*" (w„) C rf;^_^(t) and 

Wn^w weakly in GSBV^{VL; M?). 
We claim that there exists km ^ +oo such that 

(7.14) vl^.rrSt)'^w weakly in G'5BV/(r2;M2). 

Then since S9^^\vl^^^{t)) C ra(i) for all m and in view of pT^ . we deduce that S'9(*)(u;) C ra(0- 
Let us prove H7.14|l . Fixed m S N, let us choose km in such a way that 

m\\w - ipk^Wi ^ 0. 

By minimality of w^'"^(t) we have for all n 

Passing to the limit in n, by lower semicontinuity we get for some C > 

l|V«L,™(Ollp + \\vL,mit)h + m\K^,mit) - ft.^ \\i<C + m\\w - ^k^ 111. 
We deduce for m ^ +cx) 

IkL.^W-v'fc^lli-o, 

which together with ||(/3fc^ ~ ""^III ^0 implies that 

^'L.m(*) w strongly in L^{n-M^). 

Since 

l|V<„,™(i)llp + \\vl^Ut)\U <C + m\\w - 111 < C + 1 
for TO large, we have that ^^{t) w weakly in GSBV^{Vl\'E?), and this proves (|7.14|l . 

Finally, let us come to point (6). Let v G GSBVP{n;R^) with S{v) C Hs, and let us fix 
ki, . . . ,ks and toi, . . . , to,- in N. By Proposition there exists u„ G yl^^ £,(ri; M^) such that 

\im£^\t){vn)^£'\t){v) 
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and 

<M«r(^^^*Ht')\ U ^«,™.))' 

where ^(a) — > 1 as a ^ 0. Since the fc^'s and the m^'s are arbitrary, we obtain that 

(7.15) limsupf^ (^^-'*H«n) \rf::,Ji)) < M(a)f' (^^('H") \ ra(i)) ■ 

Let us suppose that wf^ J^t) Ua{t) weakly in GSBVP{Q; R^) along a suitable subsequence which 
we indicate by the same symbol. By the minimality property l|7.5|) . comparing ^{t) with ?;„ 
we get 

(7.16) s^'mutJt)) < r\t)M + £^ {s^""^'\v„)\rtjt)) +o„, 

with On — !■ as n ^ +00. Then we have 

n ' 

< hm^sup [£''{t)M + (^^-(*)(««) \ rf:,Jt))) 
< £'='(t)(i;) +limsup£^ (s''"-^'Hvn)\rtAt)) 

<£^'(i)(t^) + Mar f^^(*)(t^)\r,(o), 



that is 1)7. 10|) holds. Choosing v = Ua{t), passing to the limsup in 17.16|l . and taking into account 
(|7.15|) we obtain that 

limsupf^'(i)(wf:,a(i)) < S^'muait)). 

n 

Since by (|7.10|l £'^^{t){ua{t)) is independent of the accumulation point Ua{t), we conclude that 
((7IT|) holds. □ 

Remark 7.3. Using Lemma 17.21 it is possible to construct an increasing family {t ^a{t) : t G 
[0,T]} and a subsequence of ((5„,£„)„gN such that points (a), (6) and (c) of Lemma f7.2l hold for 
every t G [0,T]. This evolution {t Fait) : i G [0,T]} can be considered as an approximate 
quasistatic evolution, in the sense that it satisfies irreversibility, but satisfies static equilibrium 
and nondissipativity up to a small error due to the fact that a is kept fixed. The presence of 
fi{a) in the minimality property l|7.10|l takes into account the anisotropy in the approximation of 
the surface energy: in fact, since a is kept fixed, the adaptive edges of the triangulations % ai^) 
cannot recover all the possible directions. The nondissipativity condition up to a small error can 
be obtained using the minimality property (|7.10l) and following TS*, Theorem 3.13] (estimate from 
below of the total energy). 

The construction of {t Tait) : t e [0,T]} is the following. If D C [0,T] is countable 
and dense, by Lemma 17.21 and using a diagonalization argument, we can find a subsequence of 
{5„, e„)„gN and an increasing family Ta{t) G TZ{Q,b\ ^Arfi), t ^ D, such that points (a), (5) and (c) 
hold for every t E D. Let us set for every i G [0, T] 

r^(t):= n Tais). 

s>t,seD 

Clearly {t r+(i) : t G [0,T]} is increasing, in the sense that Tais) CTa{t) for all s < t. As a 
consequence, the set J of discontinuity points of 7i^(r+(i)) is at most countable. We can extract 
a further subsequence of ((5„,e„)„gN such that Ta{t) is determined also for alH G J (notice that 
ra{t)cr+{t)). ForaUi G [0, r]\(i:>U J) we set ra(i) :=r+(i). We have that Fa (t) G 7^(^s; 5^17) 
and {t ^Ta{t) : t G [0, T]} is increasing. 
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For t E D U J, Ta (t) satisfies by construction points (a) , (6) and (c) of Lemma 17.21 Let us 
consider the case t G [0, T]\{DUJ). 

Concerning point (a), we have that S^^^^w) ClTa{s) for every s £ D n[t,T], so that passing to 
the intersection we get S^^*\uait)) C Ta{t). 

As for point (b), considering s e Dn [0, i[, for every v £ GSBVP{fl;M.'^) with S'(u) cUb, we 
have there exists w„ G Af^^ (j(ri;R^) such that 

n 

and 

hmsupf^ (S^'^^^K) \rf:;_,(i)) < limsupr f^^-^'^K) \ rf:;,Js)) 

<f,{a)£^ (5^(*)(«)\r,(s)). 

Then by minimahty property (|7.5|l and passing to the hmit in n we have 

£'\t){u{t)) < £''\t){v) + ^i{a)£' {s^^'\v)\Ta{s)) . 

Letting s ^ i we get that (|7.1U|) holds. Reasoning as in Lemma [7?^ we get that also (|7.11|) holds. 
Finally, coming to point (c), we have that for all s G -D fl [0, t[ 

YirniniriVilJt)) > ]imM£^Tl:Js)) > £^{ra{s)), 

n ■ n 

SO that letting s y t, and recalling that i is a continuity point for £'^(rf^ a(')); obtain that the 
lower semicontinuity holds. 

We can now let a ^ 0. 

Lemma 7.4. There exist a map {t T{t) G Tl{^lB',dN^),t G [0,T]} and sequences 6n — > 0, 
e„ — > 0, a„ such that the following facts hold: 

(a) F" C F(s) G F(t) for allO<s<t< T; 

(b) for all t G [0,T], if Wn G ^f^^_^(fi; R^) with S'^'"^'^\wn) Q T^^l^a„{t) is such that 

Wn^w weakly in GSBVP{n; R^)^ 

then we have 

Ss('\w)cT{t); 

(c) for all t G [0,T] and for every accumulation point u{t) of (u^^ a„i^))neN for the weak 
convergence in GSBVP{n;R'^) and for all v G G'5BVP(f^; R^) with S{v)cTIb, we have 

(7.17) S^'muit)) < £'^\t){v) + £^ (^^(*)(«) \ F(i)) , 
and 

(7.18) r\t){u{t))^\imr\t){uil^,Si))- 

(d) for all t G [0, T] we have 

(7.19) £:^(F(t))<liminff^(Ff::,„Ji))- 

Proof. Let us consider 5^ ^ and £h 0. Given a g]0, ■i[ and t G [0,T], let Ta{t) be the 
rectifiable set given by Lemma 17?^ Recall that by (|7.12|) we have 
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where v'^ ^{t) is the weak hmit in GSBV^{n; R^) along a suitable subsequence depending on a of 
a minimum of the problem 

(7.20) mm{\\S/v\\p+\\v\\, + m\\v-^k\\i : f G V^^}, 
where {(pk)keN C i^(ri;R^) is dense in L^{Q;'R'^) and 

Let a„ 0, and let _D {tj : j G N} C [0,T] be countable and dense with Q E D. Using a 
diagonal argument, up to a subsequence of {Sh, Sf,.)heNj we may suppose that for alH G 13 and for 
all n 

<™ (t) ^ <r™(0 weakly in GSBVP{n; R^). 
Moreover, we may assume that for alH e Z3 and for all n 

^lajt) "a„(i) weakly in GS'BT/P(f7; R^) 

with 

By Lemma we have that Ua^{t) satisfies the minimality property H7.1()|l . 

Up to a subsequence of (a„)„gN, we may suppose that for all fc, m and t E D we have 

(7.21) yt:„^{t) ^ VkMt) weakly in G5WP(r!;R2), 
and 

(7.22) ua„{t)~^u{t) weakly in G5'WP(f2;R^). 
For alH G _D, let us set 

(7.23) m:=\JS^^'HvkMt))- 

k,m 

By Proposition l4.2l in view of the minimality property 17.10|) and taking into account that fj,{an) — *■ 
1, we have that for aU v G GSBVP{n;M.^) with S{v)cTIb 

(7.24) £'^'{t){u{t))<£'^'{t){v)+£%S3^'\v)\m), 
and as a consequence, we obtain 

We now perform the following diagonal argument. Choose Shg, Sh„ in such a way that 

\\VoT"ito) - <o(io)||l + hklaoito) - uMWl 

+ \£''{to){ul:iao{to))~£''{to){Uao{to))\ < L 

Supposing to have constructed dh„, e/i„, we choose ^ft„+i , £h„+i in such a way that for all fc < n + 1, 
m < n + 1 and for all ti with l<i<n + lwe have 

ht^^'^-^^iU) - vl2^{U)\\l + \\ut:%a,^+AU) - «a„+.(t.)||l 

+ \£''mut::^ia„,Au))-£'\u){na„,Am < 

Let us set (5„ := Sh„ and e„ := eh„, and let us prove that T{t) defined in H7.23|l satisfies the 
properties of the Lemma. We have immediately that T{t) G TZ{flB;dN^)- 
Concerning point (d), notice that 

rf:,a„(i) = U ^''"**n<;f "W), m = U s^^'HvkMm, 

m^k 771, k 

and that for all k^m 



V 



k.m 



(t) ^ Vk,m{t) weakly in GSBVP{^; K'); 
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then H7.19|l is a consequence of Theorem 12 .11 In particular, by (|7.4|l . we get that 
(7.25) H\T{t))<C'. 

Let us come to point (6). Let w„ G ylf^_„(f7; M^) ^ith S'9'"(*)(w„) C rf^_a„(i) be such that 

^ w weakly in GSBV^iVL; R^). 
For every m G N, let us choose km in such a way that 

rn\\w - '^fc,„||i 0. 
By minimality of (0 we have for all n 

iiv^t:™ wiIp + ii<:',i" wii. + ^ikt:™ w - iii < ii v^'-iif + ii^'-ii^ + ^w^n - iii. 

By construction of and in view of 17.21|l . we have 

'^fci;™ (0 ^ weakly in G^BT^f (f!; R^). 

Then passing to the limit in n, by lower semicontinuity we get for some C > 

||VWfe„,TO(t)||p + \\Vhr^.„X^)\q^rn\vu^,rn[i) -V'fc^lll < C + m||w-V3fc„||i. 

We deduce for m ^ +cxd 

||«fc„,m(i) - "^fc^lll ^ 0, 

which together with ~ w^Ili ^0 implies that 

Vk^.mit) w strongly in L"'"(ri; M^). 

Since 

l|Vwfc„,™WII < C + m||w-^fc„J|i <C+1 

for m large, we have that 

vfc,„,m(t) ^ w weakly in G'5BV/(r2; M^). 

Since S'f (*H«fe„,m(t)) C r(t) for all to, and since H^Tit)) < C , we deduce th at Sa ^^Hw) cTjt). 

Coming to point (c), we have that 17.1811 holds by construction. Moreover H7.17II holds in view 
of (|7:^ and by the fact that weakly converges in GSBVP{n;M.'^) to u{t) defined in 

In order to prove point (a), notice that if s < t with s,t € D, we have for all k,m that 
and 

<;f "(s) + gi:{t) - gll{s) - + g{t) - g{s) weakly in GSBVP{n;R^), 

where ^^"-^ "^k.mis) are defined in H7.20|l and (|7.21|) . By point (6) we deduce that 

5^^*^ (s) + - .9(s)) = («fc.„(s)) C r(i). 
Then by the definition of T{s) we get r(s) cr(i). Finally, by the same argument, we deduce 

rocr(s). 

In order to deal with all t G [0, T], we proceed as in Remark l7.3l For all t <E [0,T]\D let us set 

r+(t):= fl r(.). 

s>t,seD 

Clearly r+(t) G 7?.(ri_B; ^atJI) and satisfies point (a), so that the set J of discontinuity points of 
H^{T~^{-)) is at most countable. We can then extract a further subsequence of e„, a„)„gN such 
that T{t) is determined also for alH G J \ D (notice that T(t) C T+{t)). For all t G [0, T]\{DU J) 
we set r(i) := r+(t). We have that T{t) G TZ{V,B;dN^) and that T{t) satisfies point (a). Let us 
see that r{t) satisfies also points (&), (c) and (d) also for i G [0, T]\{D Li J). 

Concerning point (6), for every accumulation point u{t) of (w^^ ^^(t))„gN for the weak conver- 
gence in GSBVP{n; R^), by the first part of the proof, we have that S'^'^*) {u{t)) C T{s) for all s G D 
with s > t, so that passing to the intersection, we get that S^^*^u{t)) (ZT{t). 
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Let US come to point (c). Let 

uj (t) (t) u{t) weakly in GSBVP{n; M^) 

along a subsequence rij ^ +cxd. Let us set Vj := r^^^. ,a„ . and :— ge^'. ■ Up to a further 
subsequence there exists Sj € Z? with Sj t, and such that setting Uj{sj) :— Ue^^. ,a„j {sj): we have 

(7.26) Wu.is,) - ui.s,)h + IS^'is.Ku.is,)) - £'^\s,)iuis,))\ ^ 0. 
We have that there exists u*{t) G GSBVP{fl;M.'^) such that up to a subsequence 

u{sj) u*{t) weakly in GSBVP{n;M.'^). 

By the minimality property H7.17|l of u{sj), for aU v G GSBVP{n;R'^) with S'(v) cHs, we have 
that 

r\s,)iuis,)) < £^\s,){v - g{t) + g{s,)) + S^S^^'Hv) \ Tiy)). 
Passing to the limit in j we have that for all v G GSBV^'ip.; K^) with S{v) C VIb 

(7.27) £'\t){u*{t)) < £^\t)iv) + £'{S3^'\v) \ T{t)). 

As a consequence of the stability of this unilateral minimality property, it follows that 

£'\s,){u{s,))^£'^'\t)(u*{t)). 

By ((T:^ we get 

Uj{sj) u*{t) weakly in GSBV^in; M?), 

and 

(7.28) £^\s,){u,{s,))-.£^\t){u*{t)). 

By H7.5|l . comparing Uj(t) with Uj{sj) — gj{sj) + gj(t), taking into account that 

^^^■(^^■)(u,(s,))cr,(s,)cr,(t), 

we obtain 

£''{t){u,it))<£^\s,){u,{s,)) + o, 
where Oj as j ^ +oo. Passing to the limit in j we have by 17.2811 

£^^{t){u{t)) < \iinm{£''\t){uj{t)) < \im sup £^\t){uj{t)) < £''\t){u* [t)). 
'■> j 

By (fTTfjl we deduce that (|7T7|) holds. Moreover we have that £^\t){u{t)) = £''\t)(u*{t)) and 
that £'^^{t){u{t)) is independent of the accumulation point u{t). Then we deduce that H7.18|l holds. 
Finally, concerning point (d), we have that for all s G I? n [0, <[ 

hminf r(rf:,,Jt)) > hminf r(rf:.,Js)) > £^(r(s)), 

SO that letting s y t we obtain 17.19|1 . The proof is now complete. □ 

We can now prove Theorem 17. II 

PROOF OF THEOREM rm Let ((5„, e„, a„)„gN and {t T{t) G TZiTlB; dN^),t G [0,T]} be 
given by Lemma r7.4l For all t G [0,T], let us set 

Unit) -.^uli^.^t), r„(i) :=rf:,„„(0- 

Let us see that it is possible to choose an accumulation point u{t) G GSBV^{Vl\ M^) of (u„(<))„gN 
such that {t (w(i),r(t)) : t G [0,T]} is a quasistatic growth of brittle fractures in the sense of 
Dal Maso-Francfort-Toader. Let us set 

Ms) (9>V(Vw„(s)),V5e„(s)) 

- T{s){un{s)) - (aj^(s)(u„(s)),g£„(s)) 

- Q{s){Un{s)) - (9a(s)(u«(s)),5e„(s))- 
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By growth conditions of W,J-,Q and by H7.4f) we have that there exists ip £ L^{0,T) such that 
"dnis) < ip{s) for aU n. Let us consider 

i?(s) := limsup??„(s). 

n 

By Theorem 5.5 and Lemma 4.11], for every s G [0,T] there exists u{s) accumulation point 
of (u„(s))„gN for the weak convergence in GSBVP{il;M.'^) such that 

i?(s) := {dW{Vu{s)),\7g{s)) 

~g{s){u{s))~{dg{s){u{s)),g{s)). 
Applying Fatou's Lemma (in the limsup version) to (|7.t)|l with s = 0, we have that 

£:(t)(w(i),r(t)) < limsup£(0)(u„(0),r„(0)) + / '&{s)ds. 

n Jo 

By Proposition we have that limsup„ f (0)(w„(0), r„(0)) = f (0)(m(0), r(0)), so that we get 

£{t){u{t),T{t)) <£{0){u{0),T{0))+ f d{s)ds. 

Jq 

Moreover, again by TT, Theorem 3.13], 

f(OKt),r(t)) >£(o)(w(o),r(o))+ fd{s)ds, 

Jo 

so that 

(7.29) £{t){u{t),r{t)) =£{0){u{0),T{Q))+ f d{s)ds. 

Ja 

We deduce that {t — * {u{t),r{t)) : t G [0,T]} is a quasistatic growth of brittle fractures: in fact 
by Lemma [7.41 we get that r(-) is increasing, and for t G [0, T] {u{t),T{t)) G AD{g{t)) and the 
static equilibrium holds; moreover the nondissipativity condition is given by (|7.29(l . 

Let us see that points (a) and (&) of Theorem 17.11 holds. By (|7.4|l . (u„(t))„gN is weakly pre- 
compact in GSBV^{Q.]1^^) for all t G [0,T]. Moreover by Lemma f7.4l everv accumulation point 
u{t) of {un{t))nm for the weak convergence in G S BV^ {Vl;^^) is such that S'f(*)(u(t)) C T{t) and 
the minimality property H7.7(l holds. Moreover we have 

r'(t)(u(i)) =limr'(t)(u„(t)). 

n 

Since £'^''{t){u{t)) is independent of the particular accumulation point u{t), we have that point (a) 
is proved. 

Let us come to point (6). Taking into account H7.18|l and H7.19|l . for all t G [0, T] we have 
E{t) < liminf £'„(t) < limsup £■„(*) < £'(0) + [ 'd{s)ds = E{t), 

" n Jo 

SO that (|7.8(l holds. Moreover we deduce that separate convergence of elastic and surface energies 
holds at any time, so that (|7.9|) is proved. The proof is now concluded. □ 

8. The strictly convex case 

In this section we assume that the function W{x,S,) is strictly convex in ^ for a.c. x (z 
and that the function F{t,x, z) is strictly convex in z for all t G [0,r] and for a.e. x G fit as a 
consequence, the elastic energy £'^^{t,v) is strictly convex in v for all t G [0, T], and a stronger 
approximation result is available. 
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Theorem 8.1. Lei 5 € M^i^i ([0, T], E^) n L«(r2; R^)) and let 

ge G W^'\[0,T],W^-Pin;R^) n L\n;R^)), g^it) G AFe{^;^^) for all t G [0,T] 

be such that for e ^ 

ge-^ g strongly in W^^'^{[0,T],W^-P{n;M.^) (1 L'^{n;M.^)). 

Let r° e r(il) be an initial crack and let r° „ be its approximation in the sense of Provosition TE~l\ 
Let us suppose that 

(8.1) W{x^-) is strictly convex for a.e. x e VL, 

F{t,x,-) is strictly convex for a.e. {t,x) G [0,T] x il. 

Given S > 0, e > 0, a G]0,i[, let {t ^ {ut,ait),Ti^{t)) : t e [0,T]} be the piecewise constant 
interpolation of the discrete evolution given by Proposition \6.1\ relative to the initial crack ^ 
and the boundary data g^. Then there exists a quasistatic evolution {t — > {u{t),T{t)) : t G [0, T]} 
relative to the initial crack F*^ and the boundary data g in the sense of Theorem \3 . IV and sequences 
Sn —>■ 0, En 0, a„ — > 0, such that setting 

Mt) ■=<:^ajt), F„(i) ■.^Tii^.jt), 

for all t £ [0, r] the following facts hold: 

(a) yu„{t) \7u{t) strongly in LP{n;M^^'^) and u„(i) u{t) strongly in L'i{n;R^); 

(b) £{t){un{t),Vn{ty) — £{t)(u{t),V{t)), and in particular elastic and surface energies con- 
verge separately, that is 

£''\t){Un{t)) ^ S-^'muit)), £'{Vn{t)) Spirit)). 

Proof. Let us consider the sequence (5„, e„, a„)„£N and the quasistatic growth of brittle fractures 
{t — ^ {u{t),r{t)) : t G [0, r]} given in Theorem 17. II Under assumptions (|8.1|l . we have that u{t) 
is uniquely determined, because by H7.7|l u(t) minimizes 

mm{£''\t){v) : v G GSBVPin;^^), S'f(*)(„) c F(t)}, 

and £'^\t){-) is strictly convex. We conclude by point (a) of Theorcm l7. ll that ^ u{t) weakly 

in GSBVP{^1;M.'^). Point (&) is a direct consequence of Theorem 17. II By the convergence of the 
elastic energy, we deduce that 

lim / W{x,Vunit))dx = / W{x,Vu{t)) dx, 
" Jn Jn 

lim / F(t,x,Un{t)) dx ^ / F{t,x,u{t)) dx. 
" Jn Jn 

By the assumption on the strict convexity of W and F we deduce by [7] 

Vu„(t) Vu{t) strongly in 2.^(17; M^^^), 

and 

Un{t) u{t) Strongly in L«(f2;R2). 
Point (a) is now proved, and the proof is concluded. □ 
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